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Abstract. We consider the time-dependent statistical distributions of diffu- 
sive processes in relaxation to a stationary state for simple, two dimensional 
chaotic models based upon random walks on a line. We show that the cumu- 
lative functions of the hydrodynamic modes of diffusion form fractal curves in 
the complex plane, with a Hausdorff dimension larger than one. In the limit 
of vanishing wavenumber, we derive a simple expression of the diffusion coeffi- 
cient in terms of this Hausdorff dimension and the positive Lyapunov exponent 
of the chaotic model. 



1. Introduction 

The statistical mechanics of nonequilibrium processes has been the subject of 
renewed attention. Computer studies have shown that the classical motions of sys- 
tems with large numbers of particles are typically chaotic || [ll], Q |2^, [3f| 
|37| , |46| , |47| , ftc| . These observations have motivated the theoretical study of low- 
dimensional model systems where one can obtain quantitative information about 
nonequilibrium processes in terms of the chaotic properties of the microscopic dy- 
namics of the systems. 

In this context, several relationships have been established between transport 
properties such as diffusion or viscosity and characteristic quantities of chaos such 
as Lyapunov exponents, Kolmogorov- Sinai (KS) entropy per unit time, and fractal 
dimensions Jl], |J, |[ Uli |H' 19 - These relationships have been obtained for 



two different types of dynamical systems: (a) deterministic systems with Gaussian 
thermostats and, (b) open Hamiltonian systems with absorbing boundaries. 

In the deterministic systems with Gaussian thermostats ^9) , the forces ruling 
the motion of the particles are modified in order to mimic a heat pump removing the 
excess energy introduced by the external fields or other nonequilibrium constraints. 
The modified dynamical system is still time-reversal symmetric but volumes are no 
longer preserved in phase space so that trajectories, in general, converge toward 
a fractal attractor which has an invariant Sinai-Ruelle-Bowen (SRB) measure, ab- 
solutely continuous with respect to Lebesgue measure in unstable directions and 
fractal in stable directions. On average, the phase-space volumes contract at a rate 
which increases with the external fields and, thus, with the dissipation of heat. 
Since the overall contraction rate is equal to the sum of all the Lyapunov expo- 
nents, the dissipation turns out to be related to the characteristic quantities of 
chaos. In this way, formulae have been obtained for the viscosity of an externally 
sheared fluid Jll| and also for the diffusion of a Gaussian-thermostated Lorentz 
gas 0. In this last case, the formula relating the diffusion coefficient to the sum 
of the Lyapunov exponents was also rigorously proved jj, D . In the thermostated 
Lorentz gas, the kinetic energy is constant (instead of the total energy) and the 
phase space of the flow has thus the dimension three. The Lyapunov exponents of 
this system are A + > > A~ with |A~| > A + for a nonvanishing field force F. In 
the three-dimensional phase space, the fractal attractor Af has a KS entropy equal 
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to the positive Lyapunov exponent /iks = A + according to Pesin's equality, and its 
information dimension is given by Young's formula: D\ = 2 + /iks/|A _ |- I n this ex- 
pression for the fractal dimension, the 2 stands for the integer dimensions of the flow 
direction and of the unstable direction which both are equal to one for the attractor 
of a flow. However, the attractor is fractal in the stable direction which contributes 
its partial information dimension d\ — /iks/|A _ | to the total information dimension 
of the attractor A F - The diffusion coefficient of the Gaussian-thermostated Lorentz 
gas at temperature T has been alternatively expressed as 

V = Hm (^fj [-\-{A F )-\+{A F )] 
= lim^) 2 [\X-(A F )\-h KS (A F )j 

= lim (^T) \ + (A f ) Ci (Af), (1) 

where the first form is a consequence of the average contraction of phase-space 
volumes due to the Gaussian thermostatting mechanism and the relation between 
the rate of contraction and the rate of entropy production in systems with Gaussian 
thermostats, the second form follows from the first one by Pesin's equality, and the 
third one is obtained by Young formula for the partial information codimension c\ = 
1 — d\ = 1 — h-Ks/\X~\ and from the fact that lim F ^ |A~(-4f)| = lim^^o X + (A F ) @, 
[]|, ||, |). In the limit of arbitrarily small external field F, the information dimension 
of the fractal attractor approaches the phase-space dimension in a way controlled 
by the diffusion coefficient, which is the basis of this kind of formulae. 

Very similar formulae have been obtained for open Hamiltonian systems with 
absorbing boundaries |n], |l7], [l9| . Such boundaries are naturally introduced in sys- 
tems of scattering type in which particles or trajectories undergo transient collisions 
before escaping out of a phase-space region of varying geometry. The absorbing 
boundaries have the effect of driving the system out of equilibrium without a mod- 
ification of the forces which can remain Hamiltonian, i.e., time- reversal symmetric 
and also volume preserving. The escape of trajectories leads to the formation of a 
fractal repeller, described by a invariant Gibbs measure which is fractal in both the 
unstable and the stable directions. For such systems, the rate of escape is given as 
the difference between the sum of positive Lyapunov exponents and the KS entropy 
|l2|, |3(J, calculated with respect to the invariant Gibbs measure on the repeller. This 
escape-rate formula was rigorously proved for Anosov maps with rectangular holes 
H 0]. For the conservative Lorentz gas with absorbing boundaries separated by a 
distance L, the Lyapunov exponents are A + > > A~ = — A + and the diffusion 
coefficient has been obtained by the following formulae [jl7[ [l^] where the quantities 
are evaluated for the fractal repeller TZl- 
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V = lim (-J [\+(K L ) -h KS {n L )] 
-) X + (TZ L ) Ci (TZl) 



X + (K L ) c K (K L ) , (2) 

where c\ = 1 — Hks / A + denotes the partial information codimension in the stable 
or unstable directions, which are equivalent. This information codimension can be 
replaced by the partial Hausdorff codimension ch of the fractal repeller in the limit 
L —> oo Eg]. According to Eq. (0), the fractal repeller fills the phase space in the 
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limit where the absorbing boundaries are sent to infinity. The way the dimension 
of the fractal repeller approaches the phase-space dimension is controlled by the 
diffusion coefficient, and leads to the relationship between the transport coefficient 
and the codimension. The formula (0) has been generalized to all the transport 



coefficients fO 



The similarities between formulae (|f|) and (|2|) have been discussed for several 
years H. In particular, a generalization of these formulae in the presence of both 



an external field and absorbing boundaries was obtained 44 . However, each of 



the two formulae depends very much on the type of invariant set, an attractor or 
a repeller, and on the invariant measure, SRB or Gibbs, in terms of which they 
are derived. Thus, since the corresponding measures and fractal structures are so 
different in the two cases, it has not yet been possible to provide satisfying reasons 
why the formulae are so similar. 

The purpose of the present paper is to introduce a unifying approach which tran- 
scends the aforementioned differences and in which a third and new formula can be 
derived which expresses the diffusion coefficient in terms of the positive Lyapunov 
exponent of the system and the Hausdorff dimension of a fractal curve that de- 
scribes a generalized hydrodynamic mode of the diffusion process. This Hausdorff 
dimension is not, therefore, the dimension of an attractor or a repeller. The nov- 
elty of the present unifying approach is that we consider an abstract fractal curve 
directly associated with the hydrodynamic modes of relaxation toward a stationary 
state such as the equilibrium state. These hydrodynamic modes exist whether the 
system is thermostated or not, on the condition that its spatial extension is large 
compared with characteristic microscopic lengths, in order to sustain a process of 
transport by diffusion. We suppose that the time evolution of the concentration 
of tracer particles can be decomposed by a spatial Fourier transform into modes 
characterized by a wavenumber k or, equivalently, a wavelength L = 2ir/k. Each 
mode describes the inhomogeneities of concentration having this specified spatial 
periodicity and the time evolution of a non-periodic concentration is obtained by 
superposition of all the modes. Each mode has an exponential relaxation rate. For 
a deterministic dynamical system, these hydrodynamic modes can be defined in 
analogy with the conditionally invariant measures |?4| by compensating the expo- 
nential decay of the mode amplitude with an appropriate renormalization at each 
time step. In the present context, the conditionally invariant measures defining 
the hydrodynamic modes are complex because of the Fourier transform. A theory 
of these hydrodynamic modes has been described elsewhere pll p3l, where these 
conditionally invariant measures are obtained as the eigendistributions of a Perron- 
Frobenius operator of the dynamical system. The main result we have to keep in 
mind here is that these complex measures are singular with respect to the Lebesgue 
measure so that they do not have density functions. Instead, they have cumulative 
functions defined by the measure of a variable set in phase space. These cumu- 
lative functions are complex, continuous and nondiffcrentiable. When plotted in 
the complex plane, they depict fractal curves of von Koch's type [|l3[ ||f], |||] with 
Hausdorff dimensions between one and two. A recent work by the present authors 
p7[ has showed that the singular character of the hydrodynamic modes accounts 
for the entropy production of irreversible thermodynamics. 

In the present paper, we shall show that the Hausdorff dimension of the frac- 
tal curves associated with the hydrodynamic modes is controlled by the diffusion 
coefficient, which leads to a new relationship between transport and chaos. We 
shall here prove our result for a whole class of multi-baker maps. This may appear 
restrictive but several important results were first proved for the multi-baker before 
being extended to more general chaotic systems. In particular, let us mention that 
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the multi-baker models have been used for the study of the nonequilibrium steady 
states and the entropy production ||, |22|, |25|, ||(| . 

The multi-baker models have the advantage of being exactly solvable so that 
they are therefore appropriate for the intensive studies of different nonequilibrium 
properties. Many of their properties are known in detail |24|, Q . 

The organization of the paper is as follows. In Sec. |2| we discuss general models 
of deterministic random walks, which include both conservative and so-called dis- 
sipative cases. In Sec. 0, the time evolution operator on the statistical ensembles 
of trajectories is introduced. The eigendistributions and their cumulative functions 
are derived. Two different types of stationary solutions to the time evolution oper- 
ator are discussed in Sec. [I| one for conservative systems, the other for dissipative 
ones. The computation of the Hausdorff dimensions of the cumulative functions of 
the eigendistributions is done in Sec. ||. Conclusions are drawn in Sec. ^. 



2. Deterministic Models of Diffusion 

The simplest examples of reversible dynamical systems with diffusive properties 
are based upon simple stochastic processes. The most commonly studied is the 
discrete symmetric dyadic random walk where independent particles are allowed to 
move on the sites of a one-dimensional discrete lattice and the time evolution is 
determined by the condition that the probability is the same for a particle to move 
to the nearest neighboring site, on the right or left |l8|, ^2|, ^3|, [4l| . 

In order to study this process as one with a phase-space dynamics, we consider 
the set of all the possible infinite (both in the past and future) trajectories of such 
particles. It consists of sequences {n-ij^-oo °f integers rii S Z with the restriction 
that the difference between the integers at successive time steps is always equal to 
plus or minus one, i.e. Vi G Z, 

rii — = ±1. 

The stochastic time evolution is thus replaced by the shift operator on those se- 
quences, 

E({n i })=K}, 

with v! i = ni + \. 

Alternatively, a convenient representation of phase space is to consider binary 
sequences of zeros and ones labeling respectively hops to the left and right. A 
trajectory of the random walker is then given by an integer coordinate representing, 
say, the position at time zero and an arbitrary binary sequence {^il^-oo' w i & 
{0, 1} Vi, coding all its (past and future) displacements. 

This phase space is a two-dimensional continuum, which is best represented by 
points on the unit square. The transposition from symbolic sequences to points of 
the unit square is given by the dyadic expansion of x and y. Let 

oo 

- = E^T< ( 3 ) 

i=0 
— oo 

V = W 

i=-l 

i.e. the point (x, y) on the unit square is coded by the bi-infinite sequence {^i\t^-oc^ 
with the positive indices (including zero) coding the x-component and the negative 
ones the y-component. 

The reason for introducing this representation of trajectories is that the time 
evolution represented by the shift operator on binary sequences is isomorphic to 
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the baker-map acting on points (x, y) of the unit square, 



f „u{ (2x,y/2), 0<*<l/2, 

^ y > \{2x-l,(y+l)/2), l/2<x<l. W 

The proof that the shift on binary sequences and the baker map are isomorphic is 
straightforward : the baker map, Eq. (||), takes (x,y) to (x',y') with the dyadic 
expansion 



Oi+l 



i=0 



V = 



E 



(6) 
(7) 



In other words, the couple (x',y') is coded by the shifted sequence {w^}^™^, with 

In order to describe the one-dimensional dyadic symmetric random walk, we need 
to add the aforementioned integer n which labels the site of the one-dimensional 
lattice where the particle is located at the current time. Accordingly, this random 
walk is isomorphic to the symmetric dyadic multi-baker map [El BIJ 



o , x (n-l,2x,y/2), < x < 1/2, 

^■■^^ \ {n + 1,2, -!,(,+ l)/2), 1/2<*<1, (8) 



acting in the phase space Z x [0, l] 2 . 

Likewise, a one-dimensional non-symmetric r-adic random walk, where a particle 
located at site n has probability pj to hop to the site n + j, with the integer 
j taking values between — (r — l)/2 and (r — l)/2 (assuming r to be oddQ) and 
P-(r-i)/2 + • • • +P(r-i)/2 = 1; can &l so be mapped on a deterministic multi-baker 
map. 

Similarly to the dyadic process, we code the trajectories of such a process by 
sequences {'^>i\t^L oal with £ {0, ... , r— 1}. Since every symbol Wi has respective 
weight Pu>i-(r-i)/2i the r-adic expansion of x and y takes a form slightly more 
complicated than Eqs. (§)-(§)■ We can write those expansions under the recursive 
form 

x{u ,UJi, . . . ) = ^ P]j-(r-l)/2] +P[ WD -( r _l)/2] x(ui,UJ 2 , ■ ■ ■), (9) 
2/(w_i,U;_2,...) = ^ P[( r -l)/2-j] + P[(r-l)/2-w_ 1 ] 2/(^-2, W_ 3 , ... )• (10) 



1 The case r even is treated similarly, with the modification that the probability po for a particle 
to remain at the same position is set to zero. 
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The r-adic multi-baker map that mimics the random process is given by 
' / r — 1 x \ 

[n — , ,P(r-l)/2V)> 

V 2 P-( r -l)/2 > 

< X < P_( r _i)/ 2) 

/ r - 3 x -p-( r -i)/ 2 



B r ■ (n,x,y) -> < 



n + 



„ , ,P(r-3)/ay+P(r-l)/2), 
* P-(r-3)/2 7 

P-(r-l)/2 < ^ < P-(r- 1)/2 +P-(r- 3)/2> 



r - 1 Z)j^_(/-i)/2Pj 

2 P(r-l)/2 



(r-l)/2 

,P-(r-l)/2 2/ + P 3 
3=-(r-3)/2 



(r-3)/2 

Pj < a; < 1. 

j=-(r-l)/2 



(11) 



The action of £?. r on points x and y with r-adic expansions given by Eqs. (P)-(|lO|) 
is easily found to yield 



B r [n, x{u) ,wx, ...), y(o;_i, w_ 2 , . • . )] = 

[n + u) -(r- l)/2, x(ut,W2, . . . ), y(u , v-i, ■•■)]• (12) 

Notice the reversed order of the p^s along the y-coordinate in Eq. (pr|). This 
ensures the time reversibility of the map : The involution T(x,y) — (1 — y, 1 — x) 
is a time-reversal operator or reversal symmetry for B r in the sense that 



ToBroT = BZ 1 . 



(13) 



The multi-baker map, Eq. (|1 1|) , is chaotic with the mean positive Lyapunov 
exponent 



(r-l)/2 



3=-(r-l)/2 



Pj In— = ZlKS > 
Pi 



(14) 



which is equal to its Kolmogorov- Sinai entropy, and the mean negative Lyapunov 
exponent 



(r-l)/2 

E Pi ln P-J < 

j=-(r-l)/2 



(15) 



both evaluated under the forward dynamics. 

The r-adic multi-baker maps ( |TT| ) constitute simplified models of the Poincare- 
Birkhoff mappings ruling the collision dynamics of a point-like particle elastically 
bouncing on hard disks fixed in the plane |^3|, ^5). In particular, the multi-baker 
models share many of the chaotic properties of these Lorentz-type billiards. 



Example 1. The dyadic biased random walk where a particle moves to the left with 
probability p_i = q and to the right with probability p\ = 1 — q is coded by binary 
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sequences which are mapped onto the unit square by the expansion 

oo i— 1 

x = u) q + ^2(JiqY[p(2u,-i), (16) 

i=l 1=0 

— oo 2 + 1 

V = w_i(l-g)+ 5^ wi(l-g) II pci-a,,,)- (17) 

i=-2 i=-l 

TTiis transposition yields the so-called reversible dissipative multi-baker map |24|, 

MM 

B* : (n, *, y) - ( ( " _ ^ ^ ~ 0<X<q, 
2 y ' \(n+l,(x-q)/(l-q),l-q + qy), q<x<l. 

Example 2. The triadic totally symmetric random walk, with p_j = po = Pi = 
1/3, was discussed by the authors in [ p7| . Similarly to Eqs. ([|)-(0j, infinite se- 
quences of zeros, ones and twos are mapped onto the unit square by the expansion 



EU>i 

4=0 



(19) 



— OO 

v = E^t- (20) 

i=-l 

The triadic multi-baker map is simply 

\n-l,3x,y/3), < x < 1/3, 

B, : (//..,-.//) - \ (n,3a;-l,(y + l)/3), 1/3 < z < 2/3, (21) 
(n + l,3x-2,(y + 2)/3), 2/3 < x < 1. 

Example 3. TTie biased triadic random walk, where a particle hops to the left with 
probability p_i = Sl, to the right with probability p\ = sr, and remains at its 
position with probability p Q — sn = 1 — «l — s Ri is one of the processes studied 
by Tel, Vollmer, and Breymann |2|, |^, ^TJ| j>he transposition of the infinite 
symbolic sequences to the points of the unit square is similar to Eqs. (f?6[)-(|7^) and 
yields the reversible dissipative triadic multi-baker map 

i(n-l,x/s L ,s R y), 0<x<s L , 
(n, (a; - s L )/s N ,s N y + s R ), s L < x < s L + s N , 

(n + l,(x — sl — sn)/sr, SLy + SN + Si?,), SL + SN < x < 1. 

(22) 



The equivalence between the non-area-preserving map j\2q ) and an area-preserving 
multi-baker map with an extra variable of energy was discussed in Ref. p3| . 

3. Time Evolution of Statistical Ensembles 

In order to study the relaxation to stationarity of a large ensemble of random 
walkers, we restrict the one-dimensional lattice to a ring of L sites and impose 
periodic boundary conditions, that is we identify site with site L. The phase 
space is thus restricted to C x [0, l] 2 with C = {1, . . . , L}, 

The time evolution of a density pt{n, x, y) is determined by the Perron- Frobenius 
operator : p t+1 (n,x,y) = Vp t {n,x,y) = p t [B~ 1 (n,x,y)]/j[B~ 1 (n,x,y)], where 
we have introduced the Jacobian of the transformation, j\B~ x {n, x, y)\ — Pi/p-i 
if Xlj>iPj — V < Tlij>iPj- The average number of particles in a cylinder Q of the 
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,th 



cell is given by jitt(f2) = // n dx dy pt(n,x,y), with time evolution /xt+i(0) 



To display the fractal forms underlying the relaxation to stationarity, we consider 
the cumulative function defined as the measure of a cylinder fi = [0, 1] x [0, y] inside 
the n th site, i. e. 



9t(n,y)= dx dy' p t (n,x,y') 



(23) 



We can assume, without loss of generality, that the distribution function pt(n, x, y) 
is uniform with respect to x, since the ir-direction is uniformly expanding. Thus 
the distribution will become uniform along the x direction, for long times, indepen- 
dently of the initial conditions. Thus we can write pt{n, [0, x] x [0, y]) = x gt{n, y) 
and find that the time evolution of the cumulative function, g t (n, y) , is given by 

/ r — 1 y 

P-(r-l)/2 9t[n + 



9t+i(n,y) 



P-(r-i)/2 9t(n-\ : 

+P-(r-3)/2 9t (n + 



2 P(r-l)/2' 

< y <P(r-l)/2) 



1 



r - 3 y ~P(r-l)/2 
2 ' P(r- 3)/2 

P(r-l)/2 < 2/ < P(r-l)/2 + P(r-3)/2, 



(r-3)/2 

XI Pj9t{n-j,l) 

j=-(r-l)/2 



+P(r-l)/2 5t 



1 J/ 2-r2=-( r -3)/2"j 
P-(r-l)/2 
(r-l)/2 

53 ft < y < L 

j=-(r-3)/2 



(24) 



In order to identify the hydrodynamic modes and their eigcnfunctions, we con- 
sider y = 1 in Eq. (p4T). This yields the matrix equation 



(r-l)/2 

ff i+ i(n,l)= 53 Pj9t(n-j,l)- 

i=-(r- 1)/2 



(25) 



The eigenmodes of this equation have the form ipk(n) = exp(ifcn) with eigenvalue 

(26) 



(r-l)/2 

Xfe = 53 Pjexp(-ifcj). 

j=-(r-l)/2 



Moreover, the values of the wavenumber fc are restricted to k = 2irm/ L, with m G Z 
(mod L), by the periodic boundary conditions. 

According to the above considerations, we may suppose that the cumulative 
function g t (n,y) of the deterministic multi-baker map may be expanded as : 



9i 



(n, v) = 53 Xfc ak ^ fc (") Fk ^ ' 



(27) 
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where ak are coefficients set by the initial conditions and Fk(y) is solution of the 
system 



< V <P( r -l)/2, 



■ exp 



' J -(r-3)/2 



exp i 



y — P{r-l)/1 
P(r-3)/2 

P(r-l)/2 < V < P(r- 1)/2 +P(r- 3)/2, 



(28) 



P(r-l)/2 
Xfc 



exp 



» ^j = -(r-3)/2» 



P-(r-l)/2 



E 



(r-l)/2 



Pj <y < 1, 



-<j=-(f-3)/2 

where we used the identity J 7 ^ ( 1 ) = 1 . 
Remark 1. As long as 

p 3 / X k < 1, Vj, (29) 

-Eg. ( pq j is a contracting functional equation of de Rham type jl^, with a unique 
continuous solution corresponding to the cumulative function of the hydrodynamic 
eigendistribution of wavenumber k of the Perron- Frobenius operator. 

On large spatial scales, the transport process in the multi-baker chain can be 
modeled by the advection-diffusion equation pa] 



d t c= -v d n c + Vd 2 n c, 



(30) 



where c = c(n, t) = <?t(n, y = 1) is the particle concentration given by the cumula- 
tive function at y = 1, v is the average velocity of the particles drifting in a biased 
random walk, and T> is the diffusion coefficient. The result (^) is a consequence 
of the decomposition (pTj) of the cumulative function in terms of the hydrodynamic 
modes ipk(n) = exp(«/crt). Indeed, if we substitute the decomposition ( p7f ) into Eq. 
(|30|), we find that the hydrodynamic modes c(n,t) — XfcV'fcM are approximate so- 
lutions under the approximation that the terms in fc 3 , fc 4 ,... are neglected, in which 
case we find that 



lnxfc = -i v k - V k 2 + 0(k 3 ) 



(31) 



This approximation is justified for the advection-diffusion equation ( J3C| ) because 
the large-scale limit is equivalent to the small wavenumber limit k — > 0. In this 
way, an expansion of the eigenvalue ( p6| ) in powers of the wavenumber k allows us 
to identify the drift velocity: 



(r-l)/2 

3 Pi 

3=-(r— 1)/2 



and the diffusion coefficient 
1 




(32) 



(33) 



in terms of the parameters of the mapping (pi]). 
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Remark 2. For a symmetric random walk, with r odd, and where pj = p-j Vj, the 
eigenvalue is real and Eq. (pqj becomes : 



(r-l)/2 

Xfe=Po + 2 ^ pj cos(kj). (34) 
i=i 



In i/us symmetric case, the drift velocity vanishes, v = 0, and the random walk is 
purely diffusive with the diffusion coefficient 



(r-l)/2 

v= E j 2 w- ( 35 ) 



for totally symmetric r-adic random walks for which pj = 1/r Vj, the eigenvalue 
is given by 



sin(fcr/2) 

Xfc = ^inW^)' (36) 



and i/ie diffusion coefficient is 



r 2 - 1 

^ = "24-- (37) 



For small fc, we can expand in powers of k and obtain the equivalent of a 
gradient expansion : 



F k (y) = F (y) + ikT(y) + O(k 2 ). (38) 



The first term of this equation is a real function Fq and corresponds to the stationary 
state of eigenvalue 1. In the next section, we will differentiate between regular and 



singular stationary solutions. The second term in Eq. (38) is purely imaginary, i.e. 
T is a real function. Note that 



— = l + ikv + 0(k 2 ), (39) 

Xk 



where v is the drift velocity (|32j) of the random walker due to a possible bias. 
Substituting Eq. (Bq) into Eq. (pa) and, using Eq. (pa), we find the general 
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expression satisfied by T(y) 



P-{r-l)/2 



\ 2 J "\ P ( r -l)/2 ) \P(r-l)/2j 



P-(r-l)/2 



(*+«) 



< y < P( r -l)/2, 



-P-(r-3)/2 



r-3 
2 



J U V P(r-3)/ 



1/1 

)/2 



+T 



/ V-P(r-l)/2 \ 
\ P(r-3)/2 y 



r(y) = < 



P(r-l)/2 < y < P(r-l)/2 + P(r-3)/2, 



+P(r-l)/2 



I \ / y-Y (r - 1)/2 D- 

) 1 I „.\ E\ / " ^j=-(r-3)/2ft 



~2~ 



y (>'-l)/2 
* ^j = -(r-3)/2 

P-(r-l)/2 



P-(r-l)/2 



Vi 



V (r_1)/2 D-<«<1 
2-y=_( r -3)/2*'J - y ^ 1 - 



(40) 

In Sec. 0, we will make the connection between this equation and the functional 



equation defining a function introduced by Takagi 38 



4. Stationary State 

The stationary eigenstate of the Perron- Frobenius operator has eigenvalue 1, with 
a cumulative function Fq given by Eq. pa) for a vanishing wavenumber k = : 



P-(r-l)/2 



F °(p(r-l)/2)' 



F (y) = { P-(r-i)/a + P-(r-a)/a p (r - 3 )/2 



etc. 



< 2/ < P( r -l)/2) 
P(r-l)/2 < y < P(r-l)/2+P(r-3)/2i 



(41) 

For a symmetric process, where Pj = p-j Vj, the solution to this equation is sim- 
ply Fo(y) — y, corresponding to the uniform density of the equilibrium stationary 
state. 

Example 4. The dyadic symmetric random walk has the stationary state given by 
Fo(y) = y- Equation ( j^Z[ ) thus becomes 



T{y) 



(42) 



y+\T{2y), < y < \, 

l-y+±T(2y-l), \<y<l. 

The solution to this functional equation is the Takagi function |3^ , , which justi- 
fies the name generalized Takagi function for T in Eq. (f^7|j. The Takagi function is 
known to be a nowhere differentiable function with H aus dor ff dimension Dn(T) = 1. 
As will be seen later, this property is shared by the whole class of functions T given 
byEq. 

For a non-symmetric process, the contraction and expansion factors differ so 
that, in this case, Fq is a Lebesgue singular function Eg). 
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0.2 0.4 0.6 0.8 1 

y 



Figure 1. For the biased random walk of Example [j], Fq is a 
Lebesgue singular function, obtained as the solution of Eq. (|43| ), 
and which we display for values of q ranging from 0.3 (bottom 
curve) to 0.7 (top curve). The straight line corresponds to q = 1/2. 
Each curve was computed at 2 10 points connected together by lines. 



Example 5. The cumulative function of the stationary state of the reversible dis- 
sipative multi-baker map, Eq. of Example^ satisfies 

, ^ f^fc)' 0< V <l-q, 

We show in Fig. [i| this function evaluated for q = 0.3, 0.35, . . . , 0.7. For this map, 
the generalized Takagi function, Eq. (Uu), is found to be a solution of 



\2 9 [i-r (,/)] + (i-5)r(ic=iJ=i), l - 9 < »< l. 

This function is depicted in Fig. [| for q = 0.5 to 0.7. For values of q lower than 
0.5, we note that T has the symmetry T[q](y) = T[l — q](l — y). 



5. Hausdorff Dimension of the Hydrodynamic modes 

In this section, we investigate the fractal properties of the cumulative functions 
Ffe of the hydrodynamic modes for k ^ 0, as given by Eq. (p8|), and we compute 
their Hausdorff dimensions. 

In order to compute the solutions of Eq. ( p8[), we return to the symbolic dynamics 
and replace y by the r-adic expansion, Eq. (|10[). We consider points y with finite 
r-adic expansion y = y(u>i, . . . ,u>i) and rewrite Eq. (|2^) in terms of the symbolic 
sequences : 

F k \y{w!, . . . = V P ^ r ^^ eX p{-ik\j - (r - l)/2]} 

+ P ±iz^m exp{-ik["i-(r-l)/2}}F k [y(<j 2 ,... (45) 
Xk 
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y 



Figure 2 . The generalized Takagi function for the biased random 
walk, Eq. ([44"|). The parameter values range from q = 0.5 which 
gives the original Takagi function, Eq. (^) (top curve at the right) 
to q = 0.7 (bottom right). Each curve was computed at 2 10 points 
connected together by lines. 



Starting with Ffe[0] = 0, we can successively compute 

Fk[v(ui)], F k [y(wi,u} 2 )], i^[y(wi,cj 2 ,w 3 )], etc. 

All these points are plotted in the complex plane where they form a fractal curve 
(Reffc , ImFk ) similar to von Koch's curve, under certain conditions (see the follow- 
ing examples). 

Example 6. Consider the triadic totally symmetric random walk of Example [|. 
From Eq. ft3£), the eigenvalues are found to be \k = (1/3) (1 + 2cosfc). The 
condition, Eq. (p^j, that Eq. ( p^ j be contracting becomes (1 + 2cosfc) > 1, which 
is satisfied as long as — ir/2 < k < ir/2. In Fig. we plot the imaginary versus 
real parts of F^ for k = 0.1 and 0.5. The first of these two graphs for the smallest 
wavenumber k = 0.1, Fig.^a, is dominated by the linear contributions to F k , i.e. 
Fj c (y) ~ y + ikT(y). As we argued in [ p7| ; the first order term in k, T, is the 
triadic equivalent of the Takagi function Eq. (j^I[). In this case, T(y) is solution of 
the functional equation 

(y+lT(3y), < y < ±, 

T(y)=U + lT(3y-l), |<2/<|, (46) 
[l-y+iT(3y-2), § < y < 1. 

The deviations - observed in FigJ^b - of the curve (Rei 7 ^, Imi 7 ^) with respect to 
the generalized Takagi function ^dj ) are due to contributions which are second order 
in the wavenumber k. As we will show, these deviations are responsible for the fact 
that the Hausdorff dimension of this curve is not one for non-zero wavenumbers k. 
In fact, D H (F k )>l. 



Example 7. Still considering the triadic totally symmetric random walk, we ob- 
serve a rather spectacular fractal for the large wavenumber k — n/3 which corre- 
sponds to the eigenvalue x-rr/3 = 2/3. This case is shown in Fig. [J. The points 0, 
1/3, 2/3 and 1 are at the vertices of a half hexagon, which is responsible for the 
apparent hexagonal symmetry of this figure. 
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Figure 4. F^/s for the triadic totally symmetric random walk. 
The curve is computed at 3 10 points, and connected together by 
lines. 



Example 8. Another interesting case corresponds to the critical value k c of k for 
which Eq. $2$) looses its contracting property, i.e. \k r = Pj in one of the conditions 
^2!\). This leads to interesting solutions of Eq. (fig J. For instance, in the triadic 
totally symmetric random walk, k c — ir/2 yields 

(-iF n/2 (3y), 0<y<l/3, 
F*/2(y) = l-i + F w/2 (3y-l), l/3<y<2/3, (47) 
{l-i + iF n/2 (3y-2), 2/3 < 1/ < 1. 

F 7r / 2 {y) is thus a non- continuous function of y. In the complex plane, the curve 
disintegrates into a disjoint set of points forming the square lattice Z 2 . Figure |^ 
displays this set constructed with the 3 10 points y{u>i, ■ ■ ■ , u>xo). Similar functions 
can be observed for the 2-adic and 5-adic totally symmetric random walks for k c — 
7r/3. In the 2-adic case, the set is a regular triangular lattice; in the 5-adic case, a 
regular hexagonal lattice. 

We turn to the computation of the Hausdorff dimension of Ffc, i.e., of the curve 
(RcFk, ImFk) drawn in the complex plane. The above examples show that the plot 
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Figure 5. F^/2 f° r the triadic totally symmetric random walk. 
The set is composed of the points of Z 2 . 



of (ReFk, ImFk) forms a continuous curve in the complex plane as long as \k\ < k c . 
However, the curve densely fills plain domains of the complex plane already for lower 
values of the wavenumber, kf < \k\ < k c . For these values, the Hausdorff dimension 
of the curve saturates at the dimension two of the complex plane so that D^(k) = 2 
for kf < \k\ < k c . At the limiting value \k\ = kf, the plot of (Rei^, ImFfc) forms a 
fractal curve similar to the Levy dragon which has a Hausdorff dimension equal to 
two || ||. 

For still lower values of the wavenumber |fc| < kf, we consider the Hausdorff 
measure, Tf(Fk), of the cylinder sets specified by the sequences {u>i, ... ,oji} : 



Wl,... ,Ul 



AF k [y(u>i, . . . ,uji)] 



(48) 



where 

AF fc [y(wi, 
with the notation 



, = F k[y(ui, ■■■ ,w + 1)] - F k [y(ui, ■ . . 



{uji, ... , U)l + 1} = 



,W;_i,u/; + 1}, u x < r - 1, 
+ 1,0}, uji = r- 1, 



(49) 



and we make the further convention that y{r — 1, . . . , r — 1, r) = 1. The quantities 
AFk[y(uJi, . . . , uji)] are the segments of the curve, (ReFfc, ImFk) at the I th step of its 
construction. In the complex plane, the I th approximant of the curve is covered by 
disks of diameter given by the absolute value of the segments AFfc[y(wi, . . . , wj)], 
appearing on the right hand side of Eq. (^) . We recall that the Hausdorff dimen- 
sion, Dft(k), of Fk is the value of d such that p2[ 



iim rtm = h d<D ^ 



(50) 
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With the help of Eq. (|45|) , it is straightforward to check that 

AF fe [y(wi, . . . ,u!i)} = AF t [^ 2 ,...,^)] 



IX* I 



Xk\ 
i 

K-(r-l)/2]- (51) 



1 ' 



i=l 



Hence, the Hausdorff measure 



i / 1 {r - 1)/2 \ 

rf(^) -—nil n^-M)/fl d = h-i3 S fi - • ( 52 ) 

|Xfe| »=1 V IXfel J = -(r-l)/2 / 

According to Eq. ( |50| ) , the Hausdorff dimension thus satisfies 

IXfcl - P_( r _l)/ 2 + +^(r-l)/2- l ^ 

Example 9. for i/ie triadic totally symmetric random walk, the Hausdorff dimen- 
sion of Fir /3j shown in Fig. 0, satisfies 

2 \DhW3) / lX D H (?r/3) 

Therefore D H (n/3) = ln(3)/ln(2). 

Remark 3. /< is not possible to solve Eq. ^5§j ) and find a general expression for 
Dn(k) for all processes. However for the case of totally symmetric r-adic random 
walks for which pj — 1/r, Vj, we can solve Eq. ( |5j[ ) and find 



( ln(r 



\k\ < k, 



D H (k) = { In(rxfc) ' 1 1 ; ' (54) 
[2, k f <\k\<k c . 

kf is the wavenumber such that \k = I/aA" where the curve starts to densely fill the 
complex plane. k c is the wavenumber such that Xk = above which the cumu- 
lative function F^ no longer forms a continuous curve in the complex plane. For 
the symmetric dyadic case r = 2, the Hausdorff dimension ( ^[ ) has been previously 
derived |39| . 

Example 10. A numerical computation of Eq. ffify for the totally symmetric tri- 
adic random walk is shown in Fig. |^. In this case, Dn(k) — ln(3)/ln(l + 2cosfc) 
for \k\ <k f = 1.1960..., while D H (k) = 2 for kf < \k\ < k c = tt/2. 

For small wavenumbers, i.e. \k\ <C kf, we can expand the absolute value of the 
eigenvalue \k , given by Eq. (|6|) , in terms of the diffusion coefficient (|3^) of the 
general r-adic random walk : 

|x fc | = l-Vk 2 + 0(k 4 ). (55) 

Notice that the absolute value of the eigenvalue does not depend on the drift velocity 
v in the linear and quadratic terms of its expansion in powers of the wavenumber 
k. 

To establish a connection between the diffusion coefficient and the Hausdorff 
dimension of the hydrodynamic modes (for which |fc| <C kf), we let 

Dn(k) = 1 + ak 2 + 0{k 4 ), (56) 

and expand 

J^pf^ = 1 + a k 2 J2pj tePj + 0(k 4 ). (57) 
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Figure 6. The Hausdorff dimension Z?h for the triadic totally 
symmetric random walk, given by Eq. ([34]) with r = 3. The dashed 
line depicts the extrapolation of D#(k) outside of the interval \k\ < 
k f = 1.1960.... 



Substituting Eqs. (|55|) and ( 



into Eq. (|5 
V 



we find 

V 
~ A+' 



(58) 



where A + is the positive Lyapunov exponent ( |l4| ) of our system. 
We can thus state the following: 

Theorem 1. For an r-adic multi-baker map with probabilities {p~( r -i)/2, 
,P( r -i)/2}> the cumulative function of the hydrodynamic mode of wavenumber 



Oik 4 ), (59) 
of this system and T> is its diffusion 



k, Eq. (28), forms in the complex plane a curve with the Hausdorff dimension 

Dn(k) = 1 + 

where A + is the positive Lyapunov exponent ([Lj 
coefficient as specified by Eq. (3S). 

As a corollary, we have the following: 

Theorem 2. Under the same conditions as in Theorem^, the diffusion coefficient 
of the system is given in terms of the positive Lyapunov exponent A + and the Haus- 
dorff dimension Dn(k) of the cumulative function of the hydrodynamic mode of 
wavenumber k according to 

D H (k) - 1 



v = y 



lim 



k 2 



(60) 



Remark 4. The expression ( |fi6[ ) bears a relation to the formula ^) giving the dif- 
fusion coefficient in terms of the Lyapunov exponent and the Hausdorff codimension 
of the fractal repeller TZl of orbits trapped in a 2-dimensional ordered Lorentz gas of 
size L [ fl9| . The analogy is best seen by writing Eq. $6(\ ) for a hydrodynamic mode 
of the smallest possible nonvanishing wavenumber k = 2tt / L in a periodic chain of 
length L, whereupon Eq. $6(\ ) takes the form 

P= H,n (|-)V[£ H (2^/L)-1], ( 61 ) 
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where Dh(2-k / L) is the dimension of the fractal curve Fk with k = 2tt/L. Although 
ak 2 in Eq. fibq) plays a role very similar to the codimension ch(1Zl) of Eq. it 
cannot be identified as such because it is the coefficient of a correction to Dh above 
1, which stems from the definition of Fk as a cumulative function. 

6. Conclusions 

In this paper, we have shown with the help of simple models of determinis- 
tic diffusion that, in the last stages of the relaxation toward the stationary state, 
the hydrodynamic modes have cumulative functions with fractal Hausdorff dimen- 
sions. For small enough wavenumbers, we were able to derive an expression for the 
Hausdorff dimensions in terms of the diffusion coefficient and positive Lyapunov 
exponent of the system. Accordingly, the diffusion coefficient can be computed in 
terms of the Hausdorff dimension and the positive Lyapunov exponent, which is 
the content of Theorem g. 

Our formula ([30|) establishes a new relationship between the diffusion coefficient 
and the characteristic quantities of chaos, which generalizes both the equation (|l|) 
from the thermostated-system approach and the equation (|J) from the escape-rate 
formalism, as anticipated in the Introduction. On the one hand, Eq. ( p0| ) remains 
applicable whether the system is thermostated (i.e., in the non-area-preserving case) 
or not (i.e., in the area-preserving case). On the other hand, no absorbing boundary 
is required in the present approach where Eq. ( |60| ) was derived. Accordingly, the 
number of particles is kept constant in systems to which Eq. (pG) applies. In 
this regard, our new formula ( |60| ) unifies the previous formulae (p and (||). The 
novelty of the present approach is to consider the Hausdorff dimension of an abstract 
curve in the complex plane, instead of the Hausdorff (or information) dimension of 
either an attractor or a repeller. This abstract curve is the cumulative function of 
the hydrodynamic mode of diffusion, which is a general concept of noncquilibrium 
statistical mechanics. 

In view of the generality of the considerations developed in the present paper, we 
conjecture that Eq. (|59| ) for the Hausdorff dimension of the hydrodynamic modes 
extends to spatially extended Axiom-A systems with two degrees of freedom which 
sustain transport by diffusion, and more generally to such chaotic systems with two 
degrees of freedom. 

Moreover, the similarity between Eq. (||) and Eq. ( |(3l"| ) of the escape-rate for- 
malism suggests a generalization of our formalism to other transport processes. As 
shown by Dorfman and Gaspard |h], |2(| , the escape-rate formalism generalizes to 
other transport phenomena by setting up a first passage problem in the space of 
the Helfand moment p8| ] . In much the same way, we can infer that the eigendistri- 
butions of the time evolution in the space of the Helfand moment have cumulative 
functions with fractal Hausdorff dimensions given by a gradient expansion analo- 
gous to Eq. (p9|), where the diffusion coefficient is replaced by the corresponding 
transport coefficient. These questions will be addressed in more detail elsewhere. 
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